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Abstract 

The extended Hubbard Hamiltonian on a bcc lattice is studied at half-filling 
and for a finite hopping between next-nearest-neighbours, in mean-field ap- 
proximation. An ionic insulating broken-symmetry phase is predicted for any 
hydrogenoid bcc solid in the density range 1.0 < r s < 2.6. The occurrence of 
an ionic phase would explain the failure to achieve hydrogen metallization at 
high pressures. Moreover, a metal-insulator transition is expected for sodium 
in the 100 GPa region. 
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Alkali metals crystallize in a bcc phase under ordinary thermodynamic conditions. In 
1935, Wigner and Huntington [|TJ proposed that even hydrogen should undergo a metal- 
insulator (MI) transition from a molecular to a monatomic bcc crystal, under high pressure 
in analogy to alkali metals. During the last years new excitement arose after several claims 
for the reach of hydrogen's metallization However, it is now out of doubt that all 

the observed phases of hydrogen are still molecular, and there is no evidence of metalliza- 
tion induced by a band overlap mechanism up to 191 GPa ||. From this point of view, 
hydrogen seems to be quite different from the halogens, for which MI transitions have been 
observed 

The nature of the high pressure A phase of hydrogen still remains unexplained, and 
several hypotheses have been advanced. The recent proposal by Baranowski || that the 
hydrogen molecules may develop electric dipole moments is challenging for different reasons: 
i) first of all, such a phase would be per se interesting, being a broken-symmetry ground state 
of the symmetric H 2 molecular system; ii) in addition, the existence of such an ionic phase 
could move a possible monatomic phase further towards higher pressures; Hi) finally, given 
the similarity of hydrogen with the alkali metals, the possible existence of an ionic instability 
should show up even in lithium or sodium under the proper thermodynamic conditions. 

Provided that such a broken-symmetry ionic crystal does exist, it would be desirable 
to fix the boundaries between the ionic and the monatomic bcc phases. Approaching the 
boundary from the bcc phase, we look for an ionic instability of the monatomic crystal, or, 
in other words, for a charge density wave (CDW) instability, commensurate with the cubic 
lattice. The existence of such a CDW has been recently observed in the ground state of 



sodium and potassium, under ordinary thermodynamic conditions [Kj. Such instabilities of 
the Fermi gas had been predicted |11|] as a consequence of the Coulomb electron-electron (e- 
e) repulsion, but do not give rise to any MI transistion, since the CDW is not commensurate 
with the lattice. Such very small effect does not prevent us from considering the alkali metals 
as simple 'free-electron' metals for most aspects. On the other hand, a nearest-neighbour 
tight-binding model on a bipartite bcc lattice gives rise to a perfectly nested cubic Fermi 



surface at half-filling, and any small e-e repulsion would drive the system towards a spin 
density wave (SDW) or towards a CDW commensurate with the lattice. 

In this Letter we show that, even without nesting of the Fermi surface, an hydrogenoid 
bcc crystal undergoes a MI transition towards a broken-symmetry commensurate CDW 
phase, for an appropriate bounded range of density values. Such a conclusion emerges from 
a careful analysis of the mean-field phase diagram for an extended Hubbard Hamiltonian, 
modified in order to take in due account the long range Coulomb interactions and the hopping 
between next-nearest-neighbours. Even for a spherical Fermi surface, the model predicts 
the occurrence of a broken-symmetry insulating ground state, provided that the nearest- 
neighbour repulsive interaction V exceeds some critical value. The latter is a function of the 
other energy scales and mainly of the on-site Hubbard repulsion U between two electrons 
sharing the same lattice site. While U is only slightly affected by any increase in density, 
V scales as a -1 , being a the cubic lattice spacing. Under high pressure, V may reach its 
critical value, giving rise to a MI transition, albeit in an intermediate density range; then, at 
very high densities, the large increase of the Fermi energy, scaling as a" 2 , would eventually 
stabilize the monatomic phase. In other words, we expect that under high pressure both 
lithium and sodium should undergo a MI transition from a simple metal to an ionic insulator. 
On the other hand, such a high density instability of the bcc crystal would suggest that, in 
order to stabilize a hydrogen monatomic phase, higher densities are required than previously 
estimated. 

The extended Hubbard Hamiltonian reads as: 

H= —ti ^2 C i(T C ja — h ^2 C ia c jT 

+U n H n il + \ X] V i:j n ia n ja/ , (1) 

i iy^j era' 

where c ia (c\ a ) denote the annihilation (creation) operators for an electron in the Wannier 
state centered on the ith site of a bcc lattice, with spin projection a E {f, |} n ia = c\ a c ia , t±, 
t 2 > 0, and (ij), ((ij)) restricting the sums over nearest and next-nearest-neighbour couples, 
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respectively. 

The Hamiltonian, Eq. ([[]), incorporates two major approximations: i) it neglects all bond- 
bond and bond-ion interactions; ii) it neglects any hopping term other than those between 
nearest or next-nearest neighbours. Approximating the Wannier states by atomic hydrogen 
ground state wave functions, for a density corresponding to r s ph 2.2, the larger bond-ion 
interaction term does not exceed the 30% of the corresponding ion-ion interaction. Of course, 
any extrapolation to higher densities would require some caution. Regarding the neglected 
hopping terms, we must notice that the insertion of such exponentially decreasing terms 
does not change the shape of the Fermi surface in a significant way. A ratio £2/^1 ~ 0.6 -i- 0.7 
allows for an almost spherical Fermi surface, up to a 3% deviation. 

The model may be solved in mean-field (MF) approximation by inserting (n^) = 5 + A ia 
and neglecting second order terms in the fluctuations 5n ia = n ia — (n ia .). Since we are looking 
for a commensurate CDW instability, we assume A icr = A cos(Q • Rj), being Q = C^-, 0, 0) 
the exact nesting vector of the Fermi surface when t 2 = 0. Working in the reciprocal lattice, 
i.e. introducing: 

4 = ^E^ k M CT , (2) 

ViV k 

with k summed over the N points inside the first Brillouin zone, and neglecting second order 
fluctuation terms, the Hamiltonian Eq. ([I]) reads (up to a constant) as: 

H MF = £ e(k)clc ka + T £ 4 +Q(T c k(T - NT A, (3) 

where T = A(U - WW), e(k) = £i(k) + e 2 (k), and: 

e,(k)=-A tl cos cos ^fj cos ^fj , (4) 

£ 2 (k)= — t 2 [cos(k x a) + cos(kyd) + cos(k z a)]. (5) 

Here W is a renormalized interaction parameter summing up all the long range Coulomb 
interactions: 

1 00 

W = -J2(-l) m+1 ZmV( m \ (6) 
8 m=l 
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where, for m = 1,2, ...oo, V^ m > denotes the interaction term for nearest-neighbours, 
next-nearest-neighbours, etc., and z m is the corresponding coordination number. 

The MF Hamiltonian, Eq. fl3|), is easily diagonalized by a canonical transformation. 
Let us introduce the spinorial notation = fcj^ , cj^Q^, with k restricted inside the 
cube \k a \ < it /a (a = x,y,z), which is exactly half the first Brillouin zone. Actually the 
transformation k — > k + Q maps such reduced zone onto the complementary half-zone. The 
Hamiltonian Eq. (Q) now reads: 

H MF = E *L^(k)*k,r - NT A, (7) 

ka 



where the 2x2 matrix h(k) is defined as: 



h(k) 



( e 1 (k)+£ 2 (k) A(U-16W)^ 



yA(U-16W) e 2 (k)-ei(k) f 
and is promptly diagonalized yielding the spectrum: 



(8) 



E ± (k)=e 2 (k)±^e 1 (ky + T*. (9) 

A gap opens between the two bands, E ± (k), whenever T > 2t 2 : in such a regime, the system 
is an insulator, and the total energy E to t is readily evaluated by summing E~ over the doubly 
occupied half-zone: 



Stat = EMk) - V^(k) 2 + T 2 ] - NT A. (10) 



kcr 



A gap equation is obtained by differentiating E tot with respect to the order parameter, A 

1 a 3 r d 3 k 1 



(11) 



16W-U 2 J (2n)^ £l{k) 2 + T 2 

A finite T always solves the latter condition at any coupling strengths but, for sake of 
consistency with the above assumption of dealing with an insulating phase, T must exceed 
the critical value, T c = 2t%. A consistent minimum for the total energy is e.g. found for 
(16W - U)/h > 3.84 if t 2 /t 1 = 0.8. 
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In the insulating phase, the Coulomb interaction is not screened by the conduction elec- 
trons, and long range contributions cannot be neglected. We should evaluate the parameters 
U, as diagonal matrix elements of the bare Coulomb interaction in the Wannier repre- 
sentation. If we assume V%j ~ 1/|R?, — Rj|, the renormalized parameter W follows from its 
definition, Eq. (g), as W ~ ^oimV, being a M the Cs-Cl Madelung constant, «m = 1.763 flI2| , 
and V = the first term in the expansion, Eq. (f|), i.e., the nearest-neighbour repulsive 
interaction. The problem is thus mapped back onto the standard extended Hubbard model 
with an effective number of nearest-neighbour sites z* = «m, to be compared with the bcc 



value z\ = 8. For the reasonable choice £2/^1 = 0.8 |13|], the U—V phase diagram for the 
CDW instability is shown in Fig. [1]. 

The boundary between the metallic and the ionic insulating phases is given by the simple 
linear relation U jt\ = —j + 2auV/ti, being 7 the minimum value of the ratio (16W — U)/ti, 
as emerging from Eq. ([11]) for T = r c = 2t2 (7 = 3.84 for £2/^1 = 0.8). In principle, an 
ionic metallic phase may exist just over the boundary, since the gap closure is due to band 
overlap, while a finite order parameter A always arises from the gap equation QTT]) . However, 
in presence of a band overlap, the total energy, Eq. fllpp, is incorrect, since the energy levels 
must be summed up to the Fermi value inside both the bands E ± . Thus, Eq. ( |TTD is not 
correct in the metallic phase and the existence of a stable broken-symmetry ground state is 
questionable in the metallic regime. Moreover, all the Coulomb interaction terms would be 
strongly screened by the conduction electrons, so that the symmetric A = ground state 
is expected to be more favoured for the metallic phase. The phase diagram is incomplete, 
since we have not taken in consideration the possible occurrence of SDW instabilities, which 
are likely to be present for U 3> V, though irrelevant for the following considerations. 

Let us first discuss the U—V phase diagram in relation with the behaviour of solid hy- 
drogen under pressure. The possible bcc phase of dense hydrogen would be a simple 'free- 
electron' metal, with an almost spherical Fermi surface and a Fermi energy comparable with 
the free electron value Ep = 1.84/r^ a.u., which also reproduces the observed Fermi energies 
for the alkali metals. In the metallic phase, neglecting any interaction term, the model Fermi 
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energy arises from the unperturbed spectrum given by Eqs. (|j) and ([5]) in terms of the param- 
eters ti, t 2 as Ep = e(kp) — s(0), where kp is the Fermi vector, akp ~ (67r 2 ) 1//3 . Comparing 
the latter with the definition for E-p in the free-electron case yields an estimate of t\ for any 
chosen ratio ^Ai, and for any fixed density. If we fix ^2/^1 = 0.8, then t\ = 0.259/r 2 a.u. 
Approximating the Wannier states by Is hydrogen wave functions, both U and V follow, 
in a.u. H , as U = §, V = — e~ 2R (]=> + y + f -R + |-R 2 )> where the nearest-neighbour 



distance R = {\/?>'n) 1 ^r s has been used. For large r s values, U/t\ ~ r 2 s) while V/t\ ~ r s , 
so that the equation of state for hydrogen in the U jt\-V jt\ phase diagram of Fig. |l| is just 
a parabola. Some possible states of dense hydrogen are reported on Fig. |] At very high 
density, the equation of state deviates from the parabolic behaviour, since V saturates for 
r s — > 0. However, the very high density limit is questionable, and must be considered as an 
extrapolation out of the range where the adopted approximations are reasonable. 

The phase diagram is not significatively altered by a change of the ratio t^/ti, since 
both the boundary line and the equation of state are shifted in the same direction and their 
relative changes compensate. 

If we rely on the emerging scenario, even in the very high density limit, then we would 
predict that hydrogen metallization requires r s < 1, since the bcc phase would be unstable 
towards an insulating ionic phase for 1.0 < r s < 2.6. The occurrence of such an ionic 
phase could explain the failure of all the attempts to reach the monatomic state at the 



currently achievable densities. In fact, as discussed by Chen et al. || and Ashcroft |15 
the occurrence of an IR active vibron mode could be justified by the presence of permanent 
dipole moments. Besides, the occurrence of any other molecular ionic phase cannot be ruled 
out by our approach, which only prevents the stability of the monatomic bcc structure for a 
given high density range. However, the high density limit r s ~ 1 is only qualitatively correct, 
as previously discussed, being the approximations out of control in this very high density 
regime. Therefore, we don't find any contrast with the Monte Carlo prediction |16[ of a 



transition towards the monatomic phase for r s w 1.3. Moreover, that numerical calculation 
evidenced the equivalence of the monatomic ground state energies in both the fee and in the 



bcc phases. Thus the ionic instability, lowering the ground state energy, should be relevant 
even if an fee structure were the most stable monatomic phase. 

While the lower bound of the ionic phase is only qualitatively determined by the present 
approach, the upper limit r s < 2.6 is much more reliable. Any hydrogenoid bcc solid 
should undergo a MI transition around that critical density value [|17|| , thus realizing an 



unexpected broken-symmetry ionic phase. At room pressure, r s = 3.94 for sodium and 
r s = 3.25 for lithium, so that according to Fig. p] both the elements are correctly predicted 
to be stable in the monatomic bcc phase. In order to reach the critical density r s rs 2.6, 
a very high pressure is required |18| - |20[| . Such a pressure could be really prohibitive for 
K, Rb and Cs. Besides, these heavier alkali elements undergo several pressure induced 
structural transitions which are believed to be driven by the electronic transfer to upper 



empty bands [18]. Therefore, our attention should be focused on the lightest alkali, since 



such transfer mechanism is negligible and the required pressure could be reached by modern 



diamond anvil cell apparatus. Lithium has been compressed up to r s rs 2.8 jnj, and a phase 
transition from bcc to fee has been observed for P = 6.9 GPa (r s rs 2.9). The occurrence of 
the fee phase could in principle invalidate our discussion, even if the ground state energies 
of such cubic monatomic structures are so close that the ionic instability cannot be ruled 
out at higher densities. Mostly, the best alkali with which we may compare our prediction 
is sodium: i) it has an almost spherical Fermi surface; ii) its 3s orbital is expected to 



be comparable for extension with the Is orbital of hydrogen ||21|| ; Hi) no structural phase 



transition has been observed up to r s rs 2.9 | 20fl ; iv) the first empty d band is far from the 



Fermi energy. A structural phase transition has been predicted |22| from bcc to hep at the 
very high density r s < 2.5, which should be reachable in the 100 GPa region fl8| . According 
to Fig. J], a MI transition should occur first, around r s rs 2.6, then the ionic phase could 
push to higher densities the structural transition. 

At this stage, we should reverse our starting question, and we should ask instead: Why 
is sodium a metal? In a broken-symmetry ionic phase, the inter-ion Coulomb interactions 
add a considerable contribution to the total ground state energy. This very same term 
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amounts to the almost entire cohesive energy of any ionic crystal. We cannot neglect such 
interactions, even at high densities, in comparison with the Fermi energy On the other 
hand, the role played by the on-site Hubbard U and by the nearest-neighbour effective 
interaction W is competitive, as it is evident from Eqs. ( |i0|) and flTTI). The former favours 
a SDW instability, whereas the latter is responsible for the onset of a CDW. In sodium, 
under normal thermodynamic conditions, such interaction terms compensate each other, 
precluding any instability The ground state is a metal, the interactions are strongly screened 
by the Fermi liquid, and give only rise to a small renormalization of the band parameters. In 
other words, the symmetric metallic phase is based on the equilibrium between competing 
interactions. Whenever we alter such an equilibrium {e.g. increasing W by compressing the 
distances) the system falls into a broken-symmetry insulating phase, where the interactions 
are no longer screened and play an essential role. Such a phase is expected for sodium 
in the 100 GPa region, for r s « 2.6. The ionic phase could explain the failure to reach 
hydrogen metallization, and would give a possible interpretation of the anomalous IR active 
vibron mode observed in the high pressure A phase. A very high density is required in 
order to reduce the role of the interactions in comparison with the kinetic electronic energy, 
and to restore symmetry. For a perfectly nested Fermi surface {t^ = 0), such a density 
is infinite, since the integral in the gap equation fllTl) diverges for T = 0. Out of nesting 
{t 2 > 0), a metallic phase arises around the U ~ 2auV region (Fig. |l|), which also explains 
the behaviour of all the alkali metals. 
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FIGURES 

FIG. 1. U-V phase diagram for a CDW instability, for £2/^1 = 0.8. The boundary between 
normal metal (M) and ionic insulator (I) is reported as a solid line. The dashed line represents the 
equation of state for a light alkali as hydrogen: the squares correspond to r s = 0.6 2.8. 
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